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Summary. Standard (Arnold-Liouville) integrable systems are intimately 
related to complex rotations. One can define a generalization of these, sharing 
many of their properties, where complex rotations are replaced by quaternionic 
ones. Actually this extension is not limited to the integrable case: one can 
define a generalization of Hamilton dynamics based on hyperKahler structures. 

1 Introduction. Standard integrable systems 

Consider an hamiltonian systems in n degrees of freedom with compact energy 
manifolds, integrable in the Arnold-Liouville sense. By definition, these can be 
isomorphically mapped to an oscillators system, i.e. (no sum on = 1, ...,n) 

Pk = -^kQk , Qk = ^kPk ; (1-1) 

hence can be mapped to a system manifestly equivariant under the standard 
torus action in R^". Obviously, as the map mentioned above is a diffeomor- 
phism, the system had to admit a torus action also in the original coordinates, 
and this is indeed guaranteed by the familiar condition of having n integrals of 
motion in involution. 

In other words, Arnold-Liouville integrable systems with n degrees of free- 
dom are intimately related to symmetry under the group U{1)® ...®U{\) = T". 

As well known, one can pass to action-angle coordinates (/, ip) , defined via 
Pk ~ \/7fe cos((/3fe), Qk — \/7fe sin((pfe); in these the evolution equations (1.1) read 
simply jfc = 0, <^fc = iJk{I), and the U{1) «) ... U{1) = [C/(l)]®" symmetry is 
again immediately apparent. 
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We can cquivalcntly consider, instead of the action-angle coordinates {I,(p), 
complex coordinates Zk = \/Tke^'^'° = Pk + iqk', each of them evolves as 

Zk = iuJkZk , (1.2) 

which of course has solution Zk{t) = e^'^''* Zk{0). 

Thus, the time evolution of an integrable system is given by a complex 
rotation with constant speed for given initial conditions - in each C"'^ subspace 
(the frequency uik depends in general on all the \zk\^ = Ik, k = 1, ...,n). 

It is natural to expect that not much would change if instead of a complex 
rotation we had a quaternionic one. We show here that indeed this is the case. 

It is more surprising that one can develop a coherent theory of dynamical 
systems which are in a way a quaternionic generalization of standard Hamil- 
tonian systems; these arc more precisely related to hyperkahler structures, and 
thus will be therefore called hyperhamiltonian. 

Hyperkahler structures are briefly introduced e.g. in [2, 4], see [3, 17] for a 
more complete discussion; they are the subject of increasing interest of mathe- 
maticians and physicists alike, see [13] for an updated overview (for the reasons 
of physical interest, see [6, 10], and [18]; for connections with integrable sys- 
tems, sec e.g. [7]). The construction of hyperhamiltonian dynamics is discussed 
in detail in [8]; see also [14] for the aspects related to quaternionic analysis [16]. 

2 Generalizations: Clifford dynamics 

Let us look at the relation between the [[/(I)]®" symmetry and integrability in 
an even more explicit way: with Zk = Pk + iqk, equation (1.2) is the familiar 

n-oscillator dynamics (1.1). With ^k "= the two-dimensional vector ^k = 
(Pk, qk) and J the standard two-dimensional symplectic matrix, this is rewritten 
as 

ik = ^k{\ii\,-,\U) J ^k ■ (2.1) 

In this notation, the integration of (1.2) goes simply through the fact that 
l^fel are constant (due to = — J) and = — /, so that exp[a;feJ] = cos{ukt) + 
isva.{ujkt). Hence the solution of (2.1), which is ^k{i) = exp[a;/; Jf]^/;(0), reads 

a(i) = [cos(a;fci)/ + sin(a;fci) J] ^(0) . (2.2) 

Once we have looked at oscillator dynamics in such an elementary way, it is 
easy to find a direct (but non entirely trivial) generalization. 

Consider indeed m-dimensional vectors S^k G (we write pk = \^k\^) and 
evolution equations of the form 

p 

ik = ^kcKai (2.3) 

a=l 

where Vka = fka{pi, Pp) are smooth functions, and Ka = Ka{pi, ■■■,pp) (with 
a = 1, ...,p) are m-dimensional matrices satisfying K]^ = —Ka and 

{Ka,Kfi} = -26a0l. (2.4) 
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We will, for later reference, denote as K the Lie algebra spanned by the K^- 
Needless to say, (2.4) states that K is a Clifford algebra, with fundamental 
quadratic form —I; see e.g. [12] or [11], chapter 11. 

It is immediate to check that under (2.3) the pk are conserved, due to K]^ = 
—Ka', thus we can consider the u^a and the Ka as constant on each trajectory 
of the system. 

The solution to (2.3) is of course £,k(t) = cxp[jykaKat]£,k{0)- When we 
evaluate the exponential, we have to take into account that, due to (2.4), 

p p 

^kaKa^kliKp = — ^ (2.5) 
a,/3=l a,/3=l 

Hence, introducing the notation 

/ p \ 1/2 p ^ 

Wfe = I ^ z^L J ; Ak := ^ Ka , 

\a=l / a=l 

SO that (2.3) is also rewritten as = Ak£,k, we have at once that 

^k{t) = [cos{ukt)I + sm{iOkt)Ak] a(0) . (2.7) 

Thus, to any matrix algebra satisfying (2.4) we can associate a generaliza- 
tion of oscillator dynamics, integrable (or, to avoid any question like "what is 
integrability?" , explicitely and elementarily solvable) by construction. We will 
also refer to these as Clifford integrable systems. 

As noted above, the ^k {t) have constant norm under Clifford dynamics; thus 
each ^k S R'" evolves on the surface of a sphere C R"* of constant radius 

Pk = |Cfe(0)|, the Pk can thus be considered as constant for any given initial data 
{^i(0), ■.■,^p(0)}, and play the role of parameters. 

We will now discuss the motion of a single vector ^k € R"*; we will assume 
Pk ^ (or the motion is trivial), and write ^, p, A, ... for ^k, Pk, Ak, for ease 
of notation. Also, as p 7^ is constant in time, we can assume p = 1 with no 
loss of generality, and the flow will take place on the unit sphere S"^~^ C R™. 

Consider ^ € S"^~^ C R™; the vectors {KiS^, KpS^} (spanning a space 
K(^)), and hence the vector belong to the tangent space TjS"""^; in general, 
however, K(^) is a proper subspace of T^S™"^ (i.e., in general p < m+1). This 
means that not all the directions of motion on S"^~^ are allowed, which is surely 
a limiting condition, and a too severe one to consider this a generalization of 
Hamilton dynamics. 

There are however some cases in which K(^) = T^S*™^^ (i.e. p = m — 1); 
these correspond to the existence of a Clifford algebra of dimension m—1 acting 
in R"*, and of course a necessary (but not sufficient) condition for this to happen 
is that the sphere S"^~^ is parallelizable; this happens only for S*^, 5*^, S*^ [11, 12]. 

Indeed, the case m = 2 corresponds to standard Hamilton dynamics (recall 
we are dealing with a single "basic block" , i.e. in this case a minimal symplectic 
subspace); in the case m = 4 we are dealing with C R'' and with the Clifford 
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(2.6) 



algebra C'£{2), well known to be isomorphic to the quaternion algebra H, and 
in the following we will concentrate on this case. The case to = 8 has not been 
explored yet as for the corresponding extension of Hamilton dynamics; note 
however that the corresponding Clifford algebra is C£{3) » H® H [11, 12], and 
does not parallelize S"^ . Thus the quaternionic case is the only full extension of 
standard Hamilton dynamics along the lines considered here. 



3 Quaternionic integrable systems 

The simplest nontrivial Clifford algebra is just the quaternion one, i.e. C£{2) = 
su{2); and here we want indeed to consider the quaternionic case. As it happens, 
this has several very special properties, shared with Ci{l) = U{1). 

As we want to set our systems in the form (2.3), i.e. in rather than in 
or in H^, we should give a representation of the quaternionic imaginary units 

k over R; this is e.g. provided by the matrices 
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(3.1) 



which satisfy - as immediate to check - the quaternion algebra 



(3.2) 



Note, for later use, that to these are associated the symplectic forms uja = 
{l/2){Ka)ijdx^ A dx^ , given explicitely by uji = dx^ A d.T^ + dx"^ A d.T^. uj2 = 
dx^ A dx^ + dx^ A dx"^, and = dx^ A dx'' + dx"^ A da;^; with f2 the standard 
volume form in R* we have Wa A = 2Q (no sum on a). 

In this case we rewrite (2.3) as 



a=l 



and the general solution is immediately recovered from (2.7). This is 

^{t) = [cos{ujt) I + sin(wi) K] ^(0) 



(3.3) 



(3.4) 



and thus, for ^(0) ^ 0, describes great circles on the sphere of radius 
Vp ~ I^C-*)!- This also means that this dynamics realizes the Hopf fibration 

Note that the are now constant matrices, and the dependence on ^ 
(actually, on p = |^p) is only through the coefficients Va- This just means, of 
course, that we have chosen a basis {Ki,K2, K^} in K w su{2). 
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Wc denote by F the algebra of smooth real functions ,f(p); this is of course 
the algebra of smooth functions / : R** ^ R which are constants on spheres S^, 
i.e. invariant under rotations in so(4) « su{2) (g) su{2). 

It is quite interesting to note that we can rewrite (3.3) in a slightly different 
form. Introduce three functions /i" : R'* ^ R, with /i" = h°'{p) G F; now 
V/i" = /"(p) and we can hence rewrite (3.3) as 

3 

i = Y,KaVh" , (3.5) 

which makes completely transparent the relation with the hamiltonian case. 

Indeed, the flow X described by (3.5) can be seen as the superposition of 
three hamiltonian flows , each of them defined by the hamiltonian h" with 
the symplectic structure associated to K^, see above. Obviously the X^ do 
not commute, but however they generate a module over F. 



4 Hyperhamiltonian dynamics 

The equation (3.5) is the starting point for the extension of this setting to the 
non-integrable case. This can be defined on an arbitrary riemannian manifold 
(M, g) of dimension 4n equipped with a hyperkahler structure. We will briefly 
recall how these are defined, and then define the associated hyperhamiltonian 
dynamics. We will work locally, i.e. on a single chart of M. 

Let (M, g) be a riemannian manidfold of (real) dimension m = 4n. Assume 
this is equipped with three complex structures {a = 1, 2, 3), i.e. three (1,1) 
tensor fields such that = —I, satisfying the quaternionic relations 

YaY0 = e„0^Y^ (4.1) 

(here and below e is the completely antisymmetric Levi-Civita tensor). 

Assume moreover that {M,g) is Kahler with respect to each of the Ya', in 
this case we say that (M, g;Ya), or M for short, is a hyperkahler manifold. 

We recall that {M,g) is Kahler with respect to Ya means that the Kahler 
form u)a is closed, dw^ = 0. The Kahler form is defined by 

uJa{v,w) := g(YaV,w) . (4.2) 

We can thus associate to each complex structure a symplectic structure Wq, 
by means of the Kahler relation (4.2); in this sense a hyperkahler structure 
(manifold) can also be seen as a "hypersymplectic" structure (manifold). 

Consider now an ordered triple of arbitrary smooth functions /i" : M ^ R; 
we associate to these a triple of vector fields by (no sum on a) 

X„ J = d/i" (4.3) 
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and define the hyperhamiltonian vector field X on M associated to the triple 
{h"} as the sum of these, 

3 

X := ^ X„ ; (4.4) 

a=l 

it is trivial to check that the Xa, and therefore X, are uniquely defined. Needless 
to say, the Xa do not commute; the relation between the Ya guarantee they 
generate a module. 

Let us now consider a local chart on M and local coordinates {x^, ...,x"*} 
on this; we write as usual di := (d/dx^). The riemannian metric g will be 
represented in coordinates by a (0,2) tensor field gij{x), the complex structure 
y by a (1,1) tensor field Y^j{x), and the symplectic form by a (0,2) tensor field 
Wij{x), i.e. uj = {l/2)Wij{x)dx^ /\dx^ . From now on we will omit to write down 
the dependence on x, for ease of notation. 

The Kahler relation (4.2) implies that Wij = gipY^^. The relation X _\co = 
dH means that Wj^X'^ = diH; as W is nondcgcncratc it has an inverse and 
we write {W'^)~^ := K. With this we can also rewrite the relation above as 

The hyperhamiltonian vector field will thus be 

X' = J2K = Y1 ^oh'^ ■ (4-5) 

a a. 

Note that we have = -g'^Yj. Indeed, K'^ = W'^ = {gY)'^, and it is 
easy to check that this is just K'^ = —Yg^^, due to Y^ = —I. Recalling that 
g^ = g, and similarly for their inverses, we get the statement. If y = —Y^ 
(which is in general not true), then K = g~^Y. 

Note also that the Ka satisfy the quaternion algebra with the natural multi- 
plication between (2,0)-type tensor fields; it is immediate to check that KagKp = 
CajS-yKj — Sa^T], with T] = g~^ the contravariant metric tensor. 

It should be stressed that hyperhamiltonian dynamics shares many proper- 
ties with standard Hamilton dynamics; among these we would like to point out 
in particular the possibility of a variational formulation [8], which is a special 
case of a general situation: Poincare vector field can be described as character- 
istic fields of a variational principle based on maximal degree forms [9] . 



5 Example. The Pauli spin equation 

The natural physical application of the extension of hamiltonian mechanics to 
the hyperkahler case concerns, of course, spin systems. 

The non-relativistic evolution equation for particles with spin one-half is 
provided by the Pauli equation. In the simplest case of a particle with spin 1 /2 
and considering only the spin degrees of freedom (we do not discuss here the 
physical meaning of this setting), this is written as 

^ = (B-S)* . (5.1) 
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Here K = Anfi/h is a dimensional constant which we set to one in the following, 
\1/ is a two-components spinor, 



* = (^!) ' V'iWeC , ||*f = 1 , (5.2) 

the real vector B is the magnetic field, with components B(f) = {Bx,By,Bz), 
and S is the vector spin operator with components the Pauli a matrices. The 
linear operator M := B • S is given by 

M = (5, ^ (5^3, 

We want now to rewrite (5.1) as an equation in R''' rather than in C^. 
In order to do so, we rewrite ip± separating their real and imaginary part as 
tp± = X±+ K±\ substituting a number by an vector, 

V'± = (g) , (5.4) 

the operator of multiplication by i is represented in R^ by the standard symplec- 
tic matrix, and we can use this to write iM as a real four- dimensional matrix 
(which we do in block notation): 



Finally, the R^ representation of the Pauli equation is given by 

dt 



where 




A 



^=M, (5.6) 



(5.7) 
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Let us define the matrices as (the R** representation of) Ki ~ ia^, 
K2 = and Ks = ias. It is immediate to check that these satisfy, by 

construction, the quatcrnionic relations. 

We can rewrite A in terms of the matrices Ka as 

A{t) = By{t)Ki + Bx{t)K2 + B,{t)K3 . (5.8) 

We see immediately from (5.7) and (5.8) that the R'' representation of the Pauli 
equation does correspond to a hyperhamiltonian system, with 

h\i,t) = {i/2)By{tMr , 

h\^,t) = {l/2)Bx{t)Ur , (5.9) 
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When B docs actually not depend on t, wc have an intcgrable system. If B 
varies with t, we have explicitely time-dependent hamiltonians the 
system is not integrable, but |^| is still constant. 

This is coherent with and actually required by the physical meaning 
of ^, which represents components of the spin wave function: its total square 
modulus = J2a must be constant, and equal to one by normalization. 

It may be interesting, in the present context, to point out that the equa- 
tion (5.1) was analyzed by Carifiena, Grabowski and Marmo in the context of 
nonlinear superposition principles [5, 15]. 

6 Some final comments 

In this note we have presented a very streamlined version of our approach. We 

would now like to briefly present some final comments concerning relation to 
the complete picture. These will be rather rough due to lack of space, see [8] 
for further detail. 

(1) . First of all, it is natural to say that a general hyperhamiltonian dynam- 
ical vector field X, see (4.3) and (4.4), on {M,g) with coordinates x^, ...,a;^" is 
integrable if there is a diffeomorphism $ : x ^ ^ (where ^ = ...,^"'") such 
that if we pass to (ps; </7f , t^l) V'l) coordinates (with s = 1, n), the flow X is 
described by (3.6) with K and h depending only on p, i.e. is a "quaternionic 
oscillator" . Note that to guarantee integrability of such a hyperhamiltonian sys- 
tem it suffices to have n independent integrals of motion (the pi, ...,/)„), rather 
than 2n as it would be the case for hamiltonian systems. 

(2) . As very briefly mentioned above, hyperhamiltonian dynamics admits a 
variational formulation, albeit a non-standard one. This is based on a double 
fibration of the enlarged phase space M := M x R (the R factor should be 
thought of as the time), i.e. M -^IB^^H with i? = B x R a smooth manifold 
of dimension 4m — 1. The variational principle is expressed by choosing an 
arbitrary domain C C B and requiring zero variation of the functional 



defined on sections $ of the bmidle M — >B under the action of vector fields 
wich are vertical for this bundle and vanish on dC. Here $* is the puUback of 
$, and d e A*"-^{M) is a form defined by d := tjj + {6ns) J^a ^'^Cc, A dt, where 
s = ±1 depending on orientation matters, the (2n — 2)-forms C,a are defined as 
= a;„ A ... A with no sum on a and, with dcja = u)a, ip ■= Y^a=i ^ Ca- 
The forms and 6 play respectively the role of the Poincare form and Poincare- 
Cartan invariant in hyperhamiltonian dynamics. 

(3). It is quite obvious that we could operate a permutation of the oja, or 
also act with a rotation on the three dimensional linear space U spanned by 




8 



these, and provided a similar operation is performed on the /i" - leave the 
hyper hamiltonian flow X invariant. 

This suggest that the triple of symplectic structures should not be seen 
as the central geometrical object, and indeed the natural invariant object char- 
acterizing the theory is the unit sphere S :^ C U. This has a counterpart, 
via the metric g and the Kahler relation, in the unit sphere in the linear space 
O spanned by the complex structures Ya', the space O is also said to be a 
quaternionic structure on M, and if the Ya are unimodular then its unit sphere 
corresponds to unimodular complex structures in {M,g) [1]. 

Note that the transformations preserving the structure ("canonical trans- 
formations") will not be those preserving the triple uja (i-e. triholomorphic 
ones), but instead those mapping the sphere S C U into itself, i.e. such that 
/* : 5 — > 5 although in general with f*{uja) ^ cua- 

This discussion shows at once - at least to the reader familiar with that 
approach - the relation of hyperhamiltonian systems, in particular integrable 
ones, with twisters theory. 
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